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Analysis  and  Synthesis  of  Spatial  Filters  That 
Have  Chebyshev  Characteristics 


1.  INTRODUCTION 

Spatial  filters,  designed  by  a new  Chebyshev  synthesis  technique  recently  de- 
veloped at  AFCRL,  are  being  made  of  stratified  dielectrics.  These  filters  have 
successfully  suppressed  wide-angle  sidelobe  structures,  particularly  in  phased- 
array  antennas  scanning  a low  sidelobe  beam  over  a limited  spatial  sector.  The 
wide-angle  sidelobes,  which  are  in  fact  grating  lobes  that  have  been  suppressed 
to  the  -12  to  -16  dB  levels,  are  reduced  to  levels  of  approximately  -20  dB  with- 
out degradation  of  the  near-broadside  radiation  characteristics  of  the  array. 

Since  the  array  face  is  flat,  the  optimal  spatial  filter  consists  of  uniform  layers 
of  dielectric  material  properly  spaced,  having  dielectric  constants  suitable  to  the 
required  3patial  discrimination. 

Such  filters  also  have  application  in  radomes,  which  usually  require  wide- 
angle  scan.  Although  the  subject  has  not  been  addressed  directly  in  this  work, 
the  equations  and  figures  given  here  can  also  be  applied  to  the  design  of  filters 
intended  for  radome  purposes. 


(Received  for  publication  13  September  1974) 
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The  principles  of  layered-dielectric  frequency-domain  filters  and  impedance 
transformers  are  now  well  established,1"5’  and  insofar  as  possible  the  techniques 
for  analysis  and  synthesis  in  this  domain  have  been  extended  to  the  spatial  domain. 
The  fundamental  difference  between  synthesis  in  the  frequency  domain  and  synt'  :- 
sis  in  the  spatial  domain  arises  because  the  transmission  coefficients  of  layers 
that  have  a high  dielectric  constant  are  strongly  frequency-dependent  but  relatively 
invariant  with  the  spatial  angle  of  incidence:  if  a wave  from  a medium  of  low  di- 
electric constant  is  incident  on  a medium  of  high  dielectric  constant,  then  for  any 
angle  of  incidence  the  wave  propagation  angle  in  the  latter  is  almost  perpendicular 
to  the  interface. 

This  distinction  suggested  a fundamental  change  in  filter  design.  The  fre- 

1-5 

quency-domain  transformers  and  filters  synthesized  by  Collin  and  others  con- 
sisted of  various  dielectric  layers  sandwiched  together.  The  spatial  domain  fil- 
ters synthesized  in  the  work  reported  here  consist  of  quarter-wave  sections  of  di- 
electric separated  by  half-wave  or  full-wave  air  spaces.  The  methods  that  Collins 

16 

outlined  in  two  basic  studies  ’ were  applied  in  their  development. 

2.  PROPAGATION  CHARACTERISTICS  OF  STRATIFIED  DIELECTRICS 
WITH  PLANE-WAVE  INCIDENCE  AT  ARBITRARY  ANGLE 

2.1  General  Formulation;  Choice  of  Coordinate  Syitems 

The  end  result  of  this  analysis  is  a technique  for  synthesizing  spatial  filters 
that  control  the  radiation  pattern  of  a phased  array  used  for  limited -sector  scan. 
For  convenience,  all  propagation  coordinates  are  referred  to  the  basic  reference 
coordinates  x.y,  z of  the  array  (Figure  1),  which  excites  a spectrum  of  plane 
waves  and  evanescent  plane  waves  as  described  in  Appendix  A.  For  purposes  of 
this  analysis,  we  need  consider  only  a single  wave  radiating  in  the  direction  of  the 
vector 

k = *kx  + yky  + 2kz  = kkQ,  (1) 

1.  Collin,  R.  E.  (1955)  Theory  and  design  of  wideband  multisection  quarter-wave 

transformers,  Proc.  IRE  43(No.  21:179-185,  Feb.  1955. 

2.  Cohn,  S.  B.  (1955)  Optimum  design  of  stepped  transmission- line  transformers, 

IRE  Tr.  MTT  MTT-3(No.  2):16-21,  Apr.  1955. 

3.  Riblett,  H.  J.  (1957)  General  synthesis  of  quarter-wave  impedance  trans- 

formers, IRE  Tr.  MTT  MTT-5(No.  l):36-43.  Jan.  1957. 

4.  Young,  L.  (1959)  Tables  for  cascaded  homogeneous  quarter- wave  trans- 

formers, IRE  Tr,  MTT  MTT-7(No.  2):233-244,  Apr.  1959. 

5.  Young,  L.  (1962)  Stepped -impedance  transformers  and  filter  prototypes, 

IRE  lr,  MTT  MTT-10(No.  5):339-359,  Sept.  1962. 

6.  Collin,  R.E.  (1960)  Field  Theory  of  Guided  Waves,  McGraw-Hill,  pp.  79-93. 
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where 

kx  = Vo  = V8ln  6o  cos  dQ), 
ky  = Vo  = ko(8infloflin  V* 

kz  = kocos  V 


2 */X. 


with  u„  and  v the  direction  cosines, 
o o 


Figure  1.  Array  Geometry  Showing 
Array  Coordinates  6,  <t>  and  Polari- 
zation Coordinates  ©,  4> 

The  aperture  field  of  the  phased  array 

is  assumed  to  be  polarized  in  the  y direction,  thus  determining  the  polarization  of 
the  incident  plane  waves.  This  radiation  is  best  described  by  the  ©,  4>,  r coordi- 
nates of  Figure  1 because  in  the  absence  of  the  dielectric  layers  the  farfield  radi- 
ation of  the  array  contains  only  an  term.  The  $ angle  is  measured  in  the  yz 
plane;  the  0 angle,  from  the  positive  x axis.  In  general,  the  vector  A1  in  ©,  *,  r 
space  is  obtained  from  the  vector  A in  x,y,  z space  by  the  dyadic  relation: 


A'  = ©Aq  + *A'#+  ?A'r 

= T • A . (2) 

where  T is  the  nine-component  dyad  written: 

T = © ( -i  sin  © + y cos  ♦ cos  0 + * sin  * cos  ©)  + 

+ *(  0 - y sin  * + z cos  #)  + 

+ ? (x  cos  ©+  y cos  * sin  0 + z sin  * sin  ©) . (3) 


This  expression,  a formalized  way  of  describing  the  relationship  between  the 
various  components  of  the  vectors  along  the  unit  vector  directions,  is  equivalent 
to  the  matrix  representation: 


[-sin  © 
0 


cos  © 


cos  ♦ cos  © 
-sin  ♦ 
cos  ♦ sin  6 


(4) 
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This  statement  of  equivalence  to  a matrix  transformation  indicates  lhat  the 
process  of  taking  the  inverse  of  a dyad  is  the  same  as  that  for  a .atrix;  and  that 
the  transformation  of  a vector  from  one  sy3tom  to  a second,  and  from  that  system 
to  a third,  is  obtained  through  matrix  multiplication  of  the  rows  and  columns  of 
the  matrix  forms  or,  equivalently,  through  scalar  multiplication  of  the  dyad  by  the 
second  transforming  dyad.  As  with  the  equivalent  matrix  operation,  dyad  multi- 
plication is  not  commutative.  Relating  the  two  syctems^n  general  by 


sin  $ 


cos  * = 


^1  - sin2 

6 cos2  0 

sin  0 cos 

0 . 

cos  0 

^1  - sin2 

0 cos2  0 

sin  0 

sin  0 

(5) 


2 2 
sin  0 cos  0 


and  specifically  restricting  the  6. 0 parameters  to  0Q , 0O  (the  direction  of  the 
main  beam  of  the  array),  yields  the  dyad  in  the  righthand  rectangular  coordinate 


system  , rQ , as  follows: 


t ■ *.  Ltt  . A-M. 

- \ } 


cos  0 


- y 


v V 

A o \ . 

T7 


+ i\ 


+ z cos  e_ 


(6) 


This  choice  of  coordinates  makes  the  unit  vector  r_  coincide  with  the  direc- 

o 


tion  of  propagation.  Thus,  rQ  = k.  Transformed  by  this  dyadic  operation  into  the 


new  coordinate  system,  the  plane  wave  (see  A^oendix  A)  is  linearly  polarized,  and 
has  the  following  electric  and  magnetic  fields: 


E'  = ^ = *oC<0o.0o)e 

,2 


-jkr0 


S'  ■ ®„V 


a / \ -jkr 

®o  C<V*o>* 


(7) 
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where  rQ  is  a distance  measured  in  t ie  direction  of  propagation,  and  C(  0Q,  0Q) 
is  a parameter  (evaluated  in  Appendix  A)  for  multiplying  each  wave. 

Throughout  this  report  the  ©Q , , r0  system  is  used  for  describing  polari- 

zation of  the  wave.  For  describing  propagation  in  the  dielectric  layers,  it  is  more 
convenient  to  resolve  the  fields  into  components  along  the  spherical  coordinate 
system  ©,  £,?.  As  before,  the  subscripted  dyad  Rq  is  usect  when  r coincides 
with  rQ . These  are  the  conventional  definitions  of  perpendicular  and  parallel  po- 
larizations for  arbitrary  incidence  of  a plane  wave  on  a dielectric  slab. 

In  general  0 , 0,  r coordinates  the  vector  A*  is  related  to  A in  the  original 
x,y,  z coordinates  by  the  dyadic  relation: 

A*  = T'  • A , (C) 

where 

T*  = 0 (x  cos  0 cos  0 + y sin  0 cos  0 - i sin  0 ) + 

+ 0 (-  x sin  0 + y cos  0 + 0 ) + 

+ r ( xu  4 y v 4 z cos  0 > . 

It  is  sometimes  convenient  to  transform  vectors  from  one  system  to  the  other. 
Unit  vectors  in  the  9.  0,  r coordinate  system  are  related  to  those  in  the  ©,  $.  r 
system  by  a rotation  -6  about  the  r,  and  So  the  vector  A*  in  the  first  system  is 
written  in  terms  of  the  vector  A'  in  the  second  system  by  means  of  the  dyadic 
notation: 

= JR  • A' . (9) 

where 

R = 0 (S  cos  5 - sin  6 ) 4 
4 0 (0  sin  v + * cos  fi  ) 4 
4 r r . 

A 

The  required  angle  6 is  obtained  by  aligning  the  unit  vector  0 parallel  to  the 
plane  of  the  dielectric  layer.  Then,  when  transformed  by  the  dyad  R • T into  the 
0,  0,  r system,  the  unit  vector  z is  perpendicular  to  0,  and 

0 • R • T • z = 0.  (10) 
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Solution  of  this  equation  yields  values  of  6 defined  by  the  follo  wing  equations: 


sin  6 = 


VT72 


cos  6 = 


-cos  6 cos  6 


Equation  (81  could  also  have  been  obtained  by  letting  the  rctatton  dyad  R 
operate  on  a vector  in  the  0,  4>,  r system,  with  the  result: 


T'  • A = R • 1 • A . 


Applying  the  transformation  in  Eq.  (8)  to  the  phased-array  fields  (see  Appen- 
dix A)  leads  to  the  following  resolution  into  two  sets  of  plane-wave  fields  as  eval- 
uated along  the  k path. 


Parallel  Polarization 


-CtWsin<*o  *jkoro 
- .-v  e ; 


B.  = <k  /“>  s 
*o  9o  ° 


Perpendicular  Polarization 


-C(eo.<t>o)cos^ocoseo  -jk0r0 


V ■ E«  <V">  > 

ro  o 


B = = 0.  <!3) 

r r 
o o 

c 

As  noted  by  Collin,  parallel-  and  perpendicular-polarized  waves  propagate 
independently  through  the  air-dielectric  slab  interface,  and  can  be  recombined  at 
the  output  of  the  layered  region  to  form  a set  of  plane  waves  in  ©Q , ♦ space,  but 
the  set  will  in  general  no  longer  be  linearly  polarized.  The  amplitude  of  the 
cross-polarized  signal  induced  by  the  dielectric  slabs  is  crucial  to  the  spatial 
filter  geometry,  and  is  easily  related  to  the  scattering  parameters  of  the  filter. 
Based  on  the  plane-wave  filter  transmission  and  reflection  coefficients,  the  plane 
wave  amplitudes  at  the  output  of  the  filter  are  given  by: 


Parallel  Polarization 

—OUT  _ p T . 
Efl  " E0O  par ' 

B°UT  " T ; 

0 Par 


Perpendicular  Polarization 


„OUT  _ p T 
% ~ \ Tprp  * 

rout  _ R T 
Be  -BooTprp' 


Sp.  — * 


where  T or  T are  the  filter  transmission  coefficients.  The  electric  fields 

par  prp 

reflected  toward  the  array  by  the  filter  interface  are  I*  times  the  incident  per- 
pendicular polarisation  and  times  the  incident  parallel  polarized  signal, 
given  by: 


Parallel  Polarization 

Perpendicular  Polarization 

p REF  _ p p 
E0  - L0opar- 

E.REF  = E.  F ; 

6 <i>Q  prp 

B REF  _ B r 
0 0o  par 

R REF  _ R r 
B0  - -B0q  prp  • 

(15) 

The  filter  parameters  IV.  T „„  and  TV.  T „ are  evaluated  in  Sec.  2.2. 
r par  par  prp  prp 

The  inverse  of  the  rotation  dyr  T R is  used  to  transform  the  scattered  field 
components  back  to  the  0,  4>,  r coordinate  system  by  changing  any  vector  A*  in 
the  8 . 0O , rQ  coordinates  to  the  vector  A1  such  that 

A'  = R'1  • A'  . (II 


where 


R~  = 0(0  cos  6 + ^ sin  6)  + 

♦(-<?  sin  6 + $ cos  51  + 


Hence. 

A = 


©/-  0 «PS  ? C°SJ-  + 0 \ + 

V 

♦ 4 /-  8 - i £°8,?  6\  + r rl  • A*  . 

V VTv  ) J 


Using  this  dyadic  transformation,  with  r taken  along  the  direction  of  propagation, 
and  using  the  definition  of  the  parallel-  and  perpendicular-polarized  incident  fields 
E0q  and  E^q,  yields  the  following  expressions  for  the  transmitted  and 
crosn-poiarized  Eqq  fields  at  the  filter  output: 


1 


( 


I 


E = 


C(«o.0o..'jkr»  [40(81nV"VC°s9°’ 


« C<V*o)e 


-3kr, 


1*0  K + *o  “cross]  . 


T - T 
par  prp 


)* 


(18) 


where  K and  K are  the  filter  transmission  coefficients  from  the  incident 

cross 

E@o-linearly  polarized  field  into  ©0-  and  cross-polarized  (4>0- polarized)  fields. 


The  reflected  fields  for  these  two  polarizations  are  given  by  the  same  expressions 

rs  T and  T „ replaced  by  the 

par  prp 

These  new  fiber  coupling  coefficients 


as  in  Eo.  (18),  but  with  the  transmission  factors  T _ and  T replaced  by  the 

par  prp 

filter  reflection  coefficients  r and  T 

par  prp 

are  defined  as  ‘J  and  $croas  • By  means  of  these  expressions  it  can  be  shown 

that  since  energy  is  conserved  for  the  two  polarizations  separately 

<T  T ’ + T r*  = l ; and  T T*  + T r* 

par  par  par  par  ’ prp  prp  prp  prp 

within  the  system  totally.  Consequently. 


1),  it  is  conserved 


“cross  “cross  "* 


* _ ♦ 

cross  cross 


= 1, 


(19) 


2.2  Wave  Matrices  and  Traiumisiinn  Throu^i  Stratified  DieJertrir  Fitters 

Quarter-wave  transformers  with  Cheoyshev  behavior  throughout  a given  fre- 
quency passband  were  synthesized  by  Collin1  in  1955.  The  process  has  since 
2-5 

been  refined.  and  there  now  exists  an  extensive  tabulation  of  transformer 
parameters  for  various  Jesign  specifications.  Spatial  filter  deBign  is  similar  in 
concept  to  bandpass  filter  design,  except  for  some  differences  that  will  be  taken 
up  when  filter  design  procedures  are  introduced. 

Collin  used  the  wave  matrix  formalism  to  derive  convenient  expressions  for 
the  transmission  properties  of  layered  impedance  sections,  and  to  describe  the 

C 

spatial  properties  of  abutting  dielectric  layers. 

The  same  formalism  will  be  used  here  to  derive 
properties  of  the  stratified  dielectric  filters. 

Figure  2 represents  two  dielectric  media, 
showing  incident  and  reflected  waves  a^  and  bj 
in  medium  1 , and  incident  and  reflected  waves 
b2  and  a £ in  medium  2 . The  incident  waves  are 
assumed  to  be  either  parallel-polarized  or  per- 
pendicular-polarized, with  no  cross -polarized 


MEDIUM  i 

(«,) 


MEDIUM  2 

( Cg ) 


Figure  2,  Interface  Between 
Two  Dielectric  Media 
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components  excited.  In  this  case  a purely  algebraic  relation  exists  between  input 
and  output.  By  the  wave  matrix,  parameters  aj  and  bj  of  medium  1 are  related 
to  parameters  ag  and  bg  of  medium  2 in  such  a way  that  cascaded  networks  or 
layers  are  conveniently  analysed.  The  input-output  parameters  rre  thus  related  by: 


(20) 


Collin1  gives  the  obvious  relations  between  the  parameters  of  this 

matrix  in  terms  of  those  of  the  conventional  scattering  matrix  (see  Appendix  B). 
and  shows  that  the  wave  matrix  for  a cascade  of  networks  is  the  product  of  the 
wave  matrices  of  each  network.  The  parameter 


(21) 


the  inverse  of  the  filter  transmission  coefficient,  is  of  particular  importance  in 
filter  design.  The  square  of  the  absolute  value  of  this  parameter,  the  power  loss 
ratio,  will  be  considered  extensively  in  the  next  section. 

The  normalizations  of  signal  parameters  required  to  convert  the  electromag- 
netic boundary  value  equations  for  parallel-  and  perpendicular-polarized  waves 
into  the  proper  format  for  wave  matrix  application  are  given  in  Appendix  B.  These 
result  in  the  following  equation  for  a junction  between  two  dielectric  materials: 


where 

r - (V~  WZ2) 

1 " (1  + z l/z2) 


and 


(22) 


and  where 
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2.  = , j 

^ k cos  6 

j 

for  parallel  polarization,  and 
cos  9 


Zj  = k, 


J 


for  perpendicular  polarization,  and  where 


zj = koVv 


sin20 


and 

S - • 

Here,  is  the  reflection  coefficient  at  terminal  1 , and  T ^ is  the  transmission 
coefficient  from  terminal  1 into  terminal  2 . The  parameter  B is  defined  for  the 
incident  medium  air.  Note  that  il  the  materials  in  sections  1 and  2 were  reversed 
in  position,  in  the  matrix  would  be  replaced  by  , and  T j would  be  replaced 
by  T2=  1-Tj. 


The  wave  matrix  for  a line  of  electrical  length  k2X  is: 


(231 


Thus,  for  a filter  section  (Figure  3)  consisting  of  a 
terminating  air-dielectric  interface  between  a dielectric 
layer  of  thickness  t , and  an  air  space  of  thickness  S , the 
wave  matrix  relating  the  parameters  . bj  below  the 
air  space  to  the  signals  at  the  output  terminal  (aj , b gl 
can  be  written: 


°2  j , ,b2  1 

vwzsjmz' 

s ' 

J 

•'I  l*> 


Figure  3.  Basic 
Filter  Section 
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which  leads  to 
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This  dielectric  layer  and  an  air  space  constitute  the  basic  element  of  the 
spatial  filters  considered  in  this  report.  The  structure  is  more  general  than 
those  considered  in  the  previous  literature  because  the  goal  here  is  to  study  spa- 
tial filtering,  not  frequency  bandpass-filtering.  It  will  be  shown  that  without  air 
spaces  it  is  difficult  to  get  any  pronounced  spatial  behavior  over  narrow  angular 
sectors. 

Equation  (25)  can  be  rewritten  in  terms  of  the  scattering  parameters  Sjj  and 
for  the  special  case  of  broadside  incidence  on  a quarter-w  ive  slab  thickness, 
as  follows: 


A 


12 


S11  e 


'V 


where 


S. 


11 


t-  1 
€+  1 


s 


12 


-2JV? 

£ + 1 * 


(26) 


The  wave  matrix  of  a filter  comprised  of  a number  of  such  sections  is  ob- 
tained simply  by  multiplying  the  filter-section  waves  in  the  sequential  order  of  the 
positions  of  their  matrices  A (S.e.tl  in  the  filter  (emitting  the  input  section  of  line 
where  appropriate).  For  example. 


N 


FILTER 


n = l 


A(V 


(?7) 


Such  generalized  filter  wave  matrices  are  used  here  in  two  ways:  w)  in  the 
Chebyshev  synthesis  procedure,  they  are  used  for  computing  the  filter  power  loss 
ratio;  (2)  with  Eq.  (B3)  they  are  used  for  deriving  the  filter  scattering  parameters 
(Tp^,.  Tpar ; and  Tp^ ) for  the  parallel  and  perpendicular  polarizations. 

These  parameters  are  then  used  in  Eq.  (18)  to  define  the  desired  and  cross- 
polarised  signals  at  the  output  of  the  filter. 
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X CHEBYShCV  FILTER  SYNTHESIS 
3.1  General  Procedure* 

Collin’s1  procedure  for  transformer  synthesis  depends  on  the  demonstrated 
properties  of  the  polynomial  expression  for  the  power  loss  ratio,  defined  by  Collin 
as  the  power  ratio  associated  with  the  inverse  of  the  filter  transmission  coefficient: 


p - AjjAjj 


Using  filter  elements  consisting  of  sections  of  transmission  line,  each  quarter- 
wavelength  long  at  the  design  frequency,  Collin  proved  that  the  power  loss  ratio  was 
an  even  polynomial  of  cos  k2t  of  degree  2n , where  n is  the  number  of  sections  in 
the  filter.  He  also  showed  that  the  power  loss  polynomial  could  be  written  as  unity 
plus  a positive  constant  times  the  square  of  Chebyshev  polynomial  T (x),  and  that 
n characteristic  impedance  values  are  sufficient  to  define  the  Chebyshev  trans- 
former. with  the  ripple  level  given  by  the  value  of  the  positive  constant. 

The  fundamental  difference  between  the  present  study  and  the  transformer 
studies  of  Collin's  and  others  is  that  the  device  studied  here  is  a spatial  filter 
rather  than  a frequency  filter.  If  the  air  spaces  Sj  were  set  equal  to  zero,  this 
analysis— like  Collin's— would  also  lead  to  power  loss  ratios  that  are  even  powers 
of  cos  kzi  . The  difference  is  that  in  the  case  of  a frequency  filter  the  kQ  = 2v f\Q 
varies  directly  with  the  frequency:  in  the  case  of  a spatial  filter,  with  t . large 
compared  with  unity.  kz^  varies  very  little  over  a wide  range.  Logically,  there- 
fore. procedures  for  designing  a spatial  filter  are  not  directly  analogous  to  those 
established  for  designing  quarter-wave  transformers. 

To  consider  the  spatial  domain— assume  dielectric  constants  of  the  filter 
media  large  enough  for  the  ^-dependence  of  the  kZj  in  the  various  layers  to  be 
neglected;  also  assume  that  the  air-dielectric  junction  characteristics  are  con- 
stant with  6.  If  the  thickness  of  each  dielectric  layer  is  a quarter- wavelength  at 
the  design  frequency  and  the  air-space  distance  S between  layers  is  fixed,  then 
the  power  loss  ratios  will  be  even  powers  of  cos(kQS  cos  0)  or  even  powers  of 
sin(kQS  cos  0).  In  such  case,  the  dielectric  slab  is  considered  a lumped  reactive 
component;  filter  synthesis  consists  in  choosing  the  magnitude  of  this  reactance. 


3.2  The  Mathematic*  of  Filter  Syntbed* 

The  first  five  Chebyshev  polynomials  are: 
Tq(x)  = 1. 

Tj(x)  = x , 
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T2(x>  = 2x2- 1 , 

T?(x)  = 4x3-3x  , 

T4(x>  = 8x4-8x2+l.  (29) 

The  optimal  properties  of  these  polynomials  are  well  known,  as  are  their  root 
locations  and  in-band  and  out-of-band  characteristics.  The  functions  oscillate 
with  amplitude  unity  throughout  the  range  and  all  have  the  value  unit  at  |xl  = 1 . 

The  polynomial  Tm(xl  has  all  of  its  m zeros  within  this  passband  range. 

An  alternative  way  of  expressing  the  general  polynomial  Tm<x)  is:: 

T <x)  = coshim  cosh  *x)  . (30) 

m 

This  expression  is  valid  for  all  x . but  is  particularly  useful  for  calculating 
the  stopband  polynomial  values.  Synthesis  of  spatial  filters  based  on  Chebyshev 
polynomials  follows  conventional  procedure,  which  begins  with  the  recognition 
that  the  power  loss  ratio  is  a polynomial  in  c\en  powers  of  the  sine  or  cosine  of 
<kQS  cos  0).  The  further  specification  that  the  polynomial  be  one  that  has  m 
double  zeros  within  the  passband  is  also  common  to  the  theory  of  filter  synthesis, 
and  has  the  result  that  the  power  loss  polynomial  can  be  set  equal  to  the  expression. 


and  £ j is  the  value  of  £ at  the  passband  edge. 

This  expression  is  unity  plus  a polynomial  of  order  2m , with  double  zeros 

2 

within  the  region  sin  £ < sin  £ . , and  with  the  maximum  ripple  A within  that  band. 
* * 

The  expression  Aj  jA^  j has  tile  minimum  value  unity  at  the  polynomial  zeros. 

3.  2. 1 TVVO-LAYER  FILTERS 

The  coefficient  A^j  of  the  wave  matrix  for  a filter  made  of  two  identical  di- 
electric slabs  of  dielectric  constant  £ and  thickness  t set  to  a quarter  wavelength 
(in  e),  computed  from  Eq.  (28),  is 

A1I  ■ jH**-5-!  •'*)•  (32> 

b21 
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and  Sjj  and  Sgj  are  computed  from  Eq.  (26).  Since  S^  is  real,  the  power  loss 
ratio  is: 


A11  A11 


,4TT  K> 

s2l! 


(33) 


To  synthesize  a two-layer  filter,  this  ratio  is  set  equal  to  the  expression 
A„AU  » (34> 


Defining  a constant 
A 


G = 


sin^ 


(35) 


then  leads  to  the  following  equation  for  the  dielectric-layer  reflection  coefficient 
-1  Wl  + G~^~ 


J11 


(36) 


Solving  for  «.  _ means  of  Eq.  (36)  and  selecting  the  quarter-wave  filter  thick- 

9 

ness  for  the  given  t , completes  the  synthesis  procedure  for  given  A and  sin  ? j . 
3.2.2  THREE-LAYER  FILTERS 

Consider  three  layers,  each  a quarter-wave  thick  at  broadside  and  arranged 
symmetrically  so  that  the  dielectric  constant  of  the  central  layer  is  £2  and  that  of 
the  two  outer  layers  is  £j  , with  all  layers  again  separated  by  the  spacing  S.  The 
A j j coefficient  of  the  resulting  wave  matrix,  viewed  at  »he  filter  input,  is 


‘11 


S21(,,S2t<2’ 


jV^je^  -Sj  j(l)  Sj  j(2)  e"j?) 
Sn(l)  e'1?(sil(2)  e&  -Sn(l)  e'i?)j  . 


(37) 


where  the  parenthetical  (1)  and  (2)  dis  inguish  the  parameters  for  the  dielectric 
layer  tj  from  those  of  layer  tg. 

After  some  manipulation,  the  power  i-'ss  ratio  can  be  put  into  the  form: 


A 1 1 A 1 1 = 1 + 


S2,(l)|4  |S2,(2)|2 


. 2 


-2Sn(l)  + 4 S,,(l)  sin  ? 


IV 


j^S,1(2)  + Sj  j(l)  Sj  j(2) 

:]*• 


(38) 
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Equating  this  to 


V sin  ?l/ 


yields  the  following  equation  for  Sjjd): 


2 r ' / \ 4cos4C  S?  (1)1/ 


| [i+sn( 


and  then 


Sn<2)  = 


2Sj j(l)  cos^?j 

[l  + sJjO) 


3.  2.  3 FOUR-LAYER  FILTERS 

The  coefficient  Aj  j of  the  wave  matrix  for  a four-layer  filter  in  which  the 
dielectric  constant  of  the  two  central  layers  is  «2  an<*  that  °*  the  two  outer  layers 
is  is: 


All  - ‘ p?(.2iC-S1,(l>Sn(2)-S?,(2)<Sn<,»Sn(2>e-2K)» 

-S11(l)e'i?(sn(2)fe2j?  -Sn(l)sJ1(2)-S11<2)+S11(!)e'2j?)  . (41) 

After  substantial  manipulation  this  becomes  the  following  expression  for  the 
power  loss  ratio: 

A..  A*  = 1+, JsS  (l)sin3?* 

11  11  |s21(i)|4  |s21(2)|4  I 11 

+ sin?  j^Sjjd)  Sn(2)  + 2S11(l)Sjl(2)  + 2Sn(2)  - 6Sn(l)j|  ,(42) 
which,  upon  equating  to 


yields  the  following  equations  for  defining  S^O)  and  j(2) 
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2Sn(l) 

' S (2)2 

JL  21 


sin3 


= 0 


Su<2>  - 


- fes2^!)  t2]^ 


2Sj  jU) + 2 


i2 


+ 488^(1)  cos2  C 


4S„«» 


(44) 


3.3  Filler  I)e*i"n 

The  equations  in  Sec.  3.2  were  derived  in  terms  of  sin  (.  Since  ? = (2rr/A)S  ccs  0, 
sin  £ is  zero  at  (S/A)  = n(0.5),  with  n = 0 excluded  as  trivial.  In  principle,  there- 
fore, there  are  a number  of  different  spacings  that  will  allow  proper  spatial  filter 
synthesis.  In  practice,  the  only  reasonable  spacings  for  most  applications  are 
S/A  = 0.5 — sometimes  S/A  - 1.0— because  larger  spacings  have  a multitude  of 
spatial  pass-  and  stopbands  that  do  not  generally  suit  the  given  requirements. 

For  S = 0.5A  , there  is  a passband  (sin  £ = 0)  centered  at  cos  0=1  (broad- 
side), and  also  at  cos  0 = 0 (endfire).  For  S = A there  are  passbands  at  broad- 
side.. at  60°,  and  at  endfire.  The  60°  passband  begins  just  beyond  40°.  Thus, 
the  basic  filter  can  be  made  so  as  to  have  synthesized  filtering  properties  from 
oroadside  to  somewhat  beyond  40°,  a spacing  that  is  appropriate  if  the  antenna 
radiation  needs  no  further  reduction  for  large  0 angles. 

F liter  synthesis  procedure  begins  with  determining  the  value  of  £ j at  the  end 
of  the  passband  and  the  value  of  some  spatial  angle  variable  £ at  which  a given 
rejection  level  is  required.  Equations  (30)  and  (31)  give  the  out-of-band  rejection 
for  Chebyshev  filters  of  the  general  type.  Figures  4 through  6 show  the  rejec- 
tion ratios  (in  decibels)  for  two-,  three-,  and  four-layer  filters  for  various 
2 

values  of  A consistent  with  rejection  ratios  of  up  to  30 dB  for  values  of 

sin  £ / sin  ? . < 10.  From  these  curves  it  is  possible  to  choose  the  value  of  the 
1 2 

passband  ripple  amplitude  A that  will  provide  a given  rejection  ratio  for  a specif- 
ic value  of  sin  £/  sin  . Since  the  minimum  passband  transmission  coefficient  is 
2 1 2 

1/(1  + A ),  the  constant  A must  be  kept  moderately  small  if  excessive  ripple  is 
to  be  avoided  in  the  passband. 

Figures  4 through  6 must  be  used  n conjunction  with  Figures  7 through  9. 

In  these  latter  figures,  the  required  reflection  coefficients  Sjj(l)  and  Sjj(2) 
[computed  for  two-,  three-,  and  four- layer  filters  in  Eqs.(36),  (40),  and  (44), 

respectively]  are  plotted  versus  the  values  of  sin  at  the  end  of  the  passband; 

2 1 
the  various  A values  that  were  used  were  chosen  to  cover  a range  that  would 

give  reasonable  ripple  \ alues  while  maintaining  good  filter  rejection. 

Some  interpretation  of  these  figures  is  necessary  to  prevent  misconceptions. 

Equation  (31)  and  Figures  4 through  6 are  for  filtei  •*  that  are  synthesised  exactly 
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like  Chebyshev  filters,  a result  not  possible  in  the  spatial  domain  because  the 
reflection  and  transmission  properties  of  each  dielectric  slab  do  vary  with  6. 
Under  this  restriction  the  filters  are  not  strictly  Chebyshev-like  in  their  be- 
havior yet  are  very  nearly  so  for  large  e . These  equations  and  curves  thus  ”ep- 
resent  the  idealized  Chebyshev  filter  that  has  constant  slab  properties,  with  val- 
ues at  all  angles  of  9 the  same  as  at  broadside.  They  aid  in  synthesis,  but  the 
properties  of  the  actual  dielectric  filter  must  be  evaluated  by  considering  the  ex- 
act wave  propagation  in  the  structure  [using  Eq.  (18)  and  the  filter  scattering  pa- 
rameters obtained  from  the  filter  wave  matrix,  Eq.  (27),  for  each  polarization]. 
Such  an  analysis  is  carried  out  in  Sec.  3.4  . Similarly,  the  values  of  the  reflection 
coefficients  plotted  in  Figures  7 through  9 are  for  an  idealized  filter  slab  with  all 
parameters  constant.  The  horizontal  lines  indicate  these  reflection  coefficients 
converted  to  the  values  of  e for  quarter- wave  slabs  at  broad;  :de. 

The  values  of  the  reflection  coefficients  and  those  of  the  dielectric  constants 

o 

of  the  filter  layers  are  functions  of  A and  sir  jjj , and  are  not  direct  functions  of 
S/X.  These  figures  are  thus  useful  for  all  S,  since  they  can  easily  show  how  the 
required  dielectric  constants  depend  on  the  choice  of  spacing. 

3.4  Remit*  Obtained  With  Two  Filter* 

The  synthesis  procedure  described  in  Sec.  3.3  has  been  used  to  derive  a num- 
ber of  two-,  three-,  and  four-element  filters.  The  contour  plots  in  Figures  10 
and  11  summarize  the  transmission  characteristics  of  several  four-element 

filters  for  transmitted  and  crossed  polarizations  in  generalized  u,v  space.  These 

2 

show  four-layer  filters  synthesized  subject  to  the  parameters  A = 0.0004  and 
sir.  = 0.15,  for  S = 0.5 X and  S = X.  Since  C is  unchanged,  this  choice  of  pa- 
rameters leads  to  filters  having  identical  layers  (£j  = 3.08  ; e2  = 15.14),  but  yield- 
ing a differen.'.  scan  performance  for  each  S value. 

The  filter  transmission  factors  whose  contour  plots  are  represented  here  are 
given  by  20  log10<K)  for  the  transmitted  polarization,  and  20  lofijo^Cro8s^  for 
the  cross  polarization,  where  the  parameters  (X)  and  (Xcro88>  are  obtained  di- 
rectly from  Eq.  (18).  The  circular  contours  around  the  origin  in  u.v  space  indi- 
"-le  in*.!  within  the  passband  at  least,  the  filter  properties  do  not  depend  strongly 
on  the  angle  0.  For  larger  6 values,  0-dependence  becomes  more  pronounced 
and  filter  behavior  is  noticeably  different,  with  somewhat  more  rejection  in  the 
H plane  (v=  0)  than  in  the  E plane  (u=  0).  The  filter  behavior  for  large  scan 
angles  (©  > 60°)  is  not  plotted  for  the  S = 0.5  X case  because  of  the  complexity  of 
the  contour  plots  and  the  generally  unsatisfactory  filter  performance  at  such 
angles. 
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Figure  10.  Transmission  Characteristics  of  Four- 
Layer  Chebyshev  Filter  [£t  = 3.08;  «2=  15.14;  A2  = 0.00 
sin  = 0. 15;  S = 0.5  X]:  (a)  Filter  Rejection  for  Trans 
mittea  Polarization  4u;~ (b)  Isolation  Between  Trans' 
mitted  Polarization  4^  and  C roared  Polarir  tion  ©, 


The  comparison  of  the  two  cases  does  illustrate  that  for  given  dielectric 
layers,  much  sharper  filter  skirts  are  possible  with  X spacing  than  with  0.5  X 
spacing.  This  can  be  an  enormous  advantage — if  the  filter  properties  beyond  about 
0 = 50°  are  not  important— because  the  availability,  machining  properties,  and 
uniformity  of  materials  of  lower  dielectric  constant  (less  than  15)  are  far  greater 
than  those  of  materials  of  higher  dielectric  constant. 

The  cross- polarized  signal  components  in  Figures  10(b)  and  11(b)  indicate 
that  these  signals  are  zero  in  both  principal  planes,  and  maximum  near  0=45°; 
and  they  are  In  all  cases  reduced  at  least  13  dB  by  the  S = 0.5  X filter,  and  18  dB 
by  the  S = X filter.  For  the  specific  application  considered  here,  this  amount  of 
filtering  is  more  than  adequate,  but  for  some  other  application  the  cross-polarized 
transmission  may  be  excessive. 

Several  filter  designs  that  have  spactngs  between  0.5  X and  X have  been  in- 
vestigated. Designing  sin?  so  that  it  varies  uctween  -sin?j  at  broadside  to 
+ sin  at  the  end  of  the  desired  radiation  sector  yields  higher  sin  ?/sin  ? j ratios 
and  thus  greater  rejection  for  given  stopband  angles.  Unfortunately,  this  choice 
of  leads  to  a requirement  for  excessively  high  dielectric  constants;  but  when 
the  A^  value  is  reduced  to  one  that  results  »n  more  reasonable  d:  electric  constants 
the  overall  filter  rejection  is  increased  only  slightly  over  that  obtained  with  the 
original  0.5  X spacing.  In  such  cases  there  is  apparently  little  to  gain  by  resorting 
to  spacings  that  are  not  0.5  X or  X . 


3.5  Results  Obtained  With  an  Array  of  ParalM-Plane  Waveguides 

The  numerical  results  of  Figures  10  and  11,  computed  by  means  of  the  basic 
formulas  derived  in  this  report,  imply  incident  fields  consisting  of  one  or  several 
uncoupled  plane  waves.  The  phased  array  is  a more  complex  distribution.  For  an 
infinitely  large  array,  the  field  is  rigorously  described  in  teims  of  a spectrum  of 
plane  waves,  some  propagating  and  some  cut  off  (Appendix  A).  If  a stratified 
filter  is  used  over  the  array  aperture,  then  the  waves  reflected  by  the  filter  are 
cross-coupled  by  the  array  face.  This  intercoupled  nature  of  the  fields  implies 
that  the  freespace  properties  of  the  filter  may  be  modified  by  its  proximity  to 
the  array.  Figure  12,  however,  shows  that  the  dielectric-layer  filter  can  be 
used  for  suppressing  sidelobcs  and  grating  lobes,  and  so  long  as  the  main  beam 
radiation  is  kept  within  the  filter  passband  there  is  relatively  little  disruption  of 
the  basic  array  properties.  It  also  shows  how  the  two  filters  of  Figures  10  and 
11  reduce  the  residual  grating  lobes  of  a special  array  that  incorporates  an  aper- 

7 

ture  control  technique  to  eliminate  the  dominant  grating  lobe.  Arrays  of  this 

7.  Mailloux,  R.  J. , and  Forbes,  G.  R.  (1973)  An  array  technique  with  grating- 
lobe  suppression  for  limited-scan  application,  IEEE  Tr.  AP  AP-21(No.  5): 
597-602,  Sept.  1973.  " 
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Figure  12.  Limited-Scan  Array  Grating  Lobe  Amplitude  Locus 

( odd -mode  scanning  no  f'lter; odd-mode  scanning 

with  four-layer  Chebyshev  filter  (S=  0.5 X); odd-mode 

scanning  with  four- layer  Chebyshev  filter  <S=  X>] 


type  cannot  eliminate  all  of  the  grating  lobes.  The  lobes  that  remain  grow  with 
array  scan  angle,  and  sometimes  get  to  the  -12dB  level  at  maximum  scan.  The 
grating-lobe  amplitude  loci  of  an  array  of  parallel-plane  waveguides  whose  inter- 
element spacings  and  aperture  width  are  both  3. ’7  X are  given  by  the  solid  curves 
ot  Figure  12.  The  maximum  scan  point  11.5*  is  denoted  by  a circle  at  each  of 
the  grating  lobe  loci.  The  broken  and  dashed  curves  are  for  the  filters  of  Figures 
10  and  1 1 , respectively. 

A comparison  of  these  curves  shows  that  even  in  the  presence  of  mutual  coup  - 
ling  of  the  elements  in  the  array,  the  rejection  ratios  for  the  filters  placed  directly 
on  the  actual  array  are  almost  identical  wi*'.i  those  of  the  isolated  filter.  The  re- 
jection properties  of  the  S = X filter  are  particularly  appropriate  to  the  given  array:, 
all  of  the  near  grating  lobes  are  reduced  by  substantial  factors.  The  remaining 
grating  lobe,  near  60®  at  maximum  scan,  appears  approximately  at  the  second 
filter  passband  and  so  is  relatively  unaltered  by  the  filter.  In  practice,  this  lobe 
is  further  reduced  by  the  projection  factor  of  the  finite  array,  and  would  appear 
at  about  the  -18dB  level  with  or  without  the  filter.  Thus,  the  filter  of  Figure  11 
reduces  all  grating  lobes  to  approximately  -18  dB,  and  reduces  all  but  the  60® 
grating  lobe  to  about  -20  dB.  This  is  accomplished  with  a four-layer  Chebyshev 
design  using  materials  of  moderate  dielectric  constant. 


v..- 


4.  CONCLUSIONS 


The  practicality  of  a new  technique  for  Chebyshev  synthesis  of  spatial  filters 
has  been  demonstrated.  The  filters  designed  to  test  the  t^-ory  have  given  good 
control  of  sidelobes,  especially  in  phased-airay  antennas,  where  they  can  be  used 
as  radomes  as  well.  They  are  particularly  valuable  in  the  limited-scan  arrays  of 
precision-approach  radar  (PAR)  systems. 
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Appendix  A 

Radiated  Field  of  Infinita  Array 

The  radiated  field  of  an  infinite  array  of  the  linearly  polarized  rectangular 
apertures  shown  in  Figure  1 is  given  in  Eq.  (Al)  in  terms  of  a single- component 
hertzian  potential.  The  main  beam  of  the  array  is  steered  to  tne  angles  0Q . d0 . 

E (?)  = - j uJ  V X [T!  (x.  y.  z)J  . 

B(r>  = v[vx  l(x.y.z)]  + k2I  (x.y.zi;  <Al) 


where 


Il(x,y,  z)  = x 11  (x.y.z). 


The  potential  function  II  is  given  by:. 


1 v V e Jkrr't‘  r rr-ve,  • 'il  )kmn'rdi'dy',(A2l 


where  D and  D are  the  interelement  spacings  in  the  x and  y directions, 
x y 


/«  - - jnn\ 

= ko\xum+yyn+  2 “k7/* 


r = xx  + yy  + zz ; 


Preceding  page  Hank 


X 


and  where  E(x,.y/)  is  the  aperture  field,  linearly  polarized  in  the  y direction:  r 
is  a position  vector  in  x,y,  z coordinates,  and  r*  is  measured  in  the  aperture 
(z'=  0). 

The  direction  cosines  u and  v for  the  specified  grating  lobe  positions  um , vn 

for  a main  beam  at  u_  = sin  cos  A and  = cos  0 sin  A are  defined  by: 
o o To  o o To  J 

um  = uo  + m/Dx' 

vn  * vo  n/Dy ' <A3' 

The  values  m and  n define  propagating  waves  for  all  integers  for  which  the 
parameter 

K = k Vl  - u*  - v*  (A4) 

ran  o m n 

is  real.  Thus,  a set  of  waves  is  defined,  some  propagating  and  some  cut  off.  The 

propagating  waves  travel  in  the  direction  of  the  allowed  k7>n  values  and  are  plane; 

for  these.  K_  reduces  to  k cos  0 The  potential  fun  non  H can  thus  be 
mn  o o x 


written  as  the  spectrum: 

n (x.y.z)  = Z <*.y.  *>. 
* m.n  mn 


and  the  field  components  of  this  plane  wave  spectrum  for  the  wave  at  um . vn  are: 


^mn 

^mn 

0. 

®ymn 

ko  umvn  ^mn* 

Eymn  = 

-uk  n , 

mn  mn 

B*mn 

= -k2  u n . 

o m k mn 

E*mn 

«k v n . 

o n mn 

(A6> 

In  conventional  phased  arrays  the  8 pacings  and  are  specified  to  exclude 

all  but  the  wave  associated  with  the  main  beam  (0Q . +Q).  but  certain  special- 
purpose  arrays  intended  for  limited  sector-scanning  allow  a small  number  of 
grating  lobes  to  propagate.  The  spatial  filters  described  in  this  report  reduce 
and  control  such  lobes. 

The  magnitude  of  the  electric  field  for  the  m,  n mode  is  given  by 


k « Vl  - u*  n . 

o m mn 

and  since  this  parameter  multiplies  each  wave  it  will  now  be  defined  as 


s.  _«  * 


Appendix  B 

Matrix  Definitions  for  riiilsi  and  rarnaratH  idir  Polarization 


Figure  2 schematically  shows  the  signal  parameters  that  are  related  by  the 
wave  matrix.  Conventionally,  the  incident  electric  field  is  chosen  as  the  refe — 
ence  signal,  with  the  wave  matrix  expressing  the  linear  relation  between  input 
and  output  parameters  as: 


<B1) 


The  coefficients  in  A are  related  to  the  conventional  scattering  matrix  co- 
efficients by  die  following  definitions: 


‘11 


‘12 


ri 


A s i 

a21  Tj  • 


A22  = T2  ■ 


r r 

1 1*  2 


(B2) 


for  the  scattering  parameters  defined: 
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(B3) 


When  the  transmission  lines  lor  ports  1 and  2 have  the  same  characteristic 
impedance,  then  T,  = ro  (reciprocity).  Otherwise,  this  relationship  does  not 

1 ^ g 

hold  unless,  as  Collin  shows,  the  signal  amplitudes  are  normalized  so  that  the 

* 

power  in  a wave  of  amplitude  a.  is  given  by  a. . Throughout  the  following  anal- 
ysis, these  amplitudes  are  defined  to  be  the  E field  parallel  to  the  d'electric  junc- 
tion, and  thus  Tj  4 T2  except  when  the  two  ports  are  located  in  layers  of  the  same  t. 

Incident  fields  for  both  polarizations  are  given  in  Eq.  (13).  Since  the  projection 
of  6 onto  the  plane  of  the  junction  involves  the  factor  V.Zq^V.0  in  the  incident  medium 
(air)  and  is  kzj/kj  in  a material  of  dielectric  constant  , where 

ki  = ko 

and 

"1  - k„V‘.  - ■*»*  "too  • <“> 


then  this  projection  factor  is  used  to  multiply  the  field  vector  that  lies  in  the  plane 
of  incidence  in  order  to  find  the  tangential  component  at  the  interface.  For  paral- 
lel polarization  the  incident  is  (kQ/w)  Eg  and  the  incident  tangential  12  is 
(kzoyk0)Eg.  For  perpendicular  polarization  the  is  the  tangential  E field  and 
the  tangential  B field  is  (-k0/w)E^. 

Similar  projections  apply  within  any  dielectric  layer.  For  parallel  polarization 
the  equations  for  continuity  of  tangential  E and  B fields  at  the  interface  between 
two  media  of  dielectric  constants  €j  and  tg,  for  incident  and  reflected  wave  ampli- 
tudes (Ej  = Eg)  and  Eg  in  medium  1 . and  Eg  and  E^  in  medium  2,  are: 

(l‘.vA,)(E,*E2).  (kZ2/k2)(E3,E4). 

*,,(*!-=»)  * ^2(es-E,)-  <BS> 
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Then,  defining  signal  amplitudes  , bj , a2  . b,,  as  (kz.  k,)  times  the  E-field 
amplitudes  through  E4  yields  the  following  equations  relating  the  tangential 
E fields  at  the  boundary: 


(al+bl)  = (a2+b2) 

and 

(al  " bl)  = (Z1  z2)(a2  ” b2>  ' (B6) 

w here 

Z1  = (k2i/kf  COS  0inc)' 

(The  cos  0jnc  is  included  to  remove  the  angle  dependence  of  for  the  trivial 
«2  = < j case  where  kz.  = kQ  cos  0inc , ) 

Following  a similar  procedure  for  perpendicular  polarization  yields  the  tan- 
gential E and  B equations 


Ej  + Eg  = Eg  + E4  , 

“2i(ei  ' £2)=  kZ2(E3  - E4).  (B71 

Then,  defining  at , bj  . a2 , and  bg  to  be  equal  to  Ej  through  E4  yields  the  same 
equation  as  for  parallel  polarization,  but  with 


cos  0. 


z.  = 
] 


inc 


(B8) 


J 


For  either  case,  Eq.  (Bl),  when  transformed  into  wave  matrix  form,  becomes. 

1 r. 


(*!)■  t‘- 


■ 1 


where 


r.  = 


(j  ~ V2) 
<’  - Vs) 


and 


1 


1 


(B9) 


T.  = 


1 + . 
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